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SUMMARY 

This paper describes the calculation of shape and scale parameters of the 
two-parameter Weibull distribution using the least-squares analysis and maximum 
likelihood methods for volume- and surface-flaw-induced fracture in ceramics 
with complete and censored samples. Detailed procedures are given for evaluat- 
ing 90-percent confidence intervals for maximum likelihood estimates of shape 
and scale parameters, the unbiased estimates of the shape parameters, and the 
Weibull mean values and corresponding standard deviations. Furthermore, this 
paper describes the necessary steps for detecting outliers and for calculating 
the Kolmogorov-Smi rnov and the Anderson-Darl i ng goodness-of-f i t statistics and 
90-percent confidence bands about the Weibull distribution. It also shows how 
to calculate the Batdorf flaw-density constants by using the Weibull distribu- 
tion statistical parameters. The techniques described have been verified with 
several example problems from the open literature, and have been coded in the 
Structural Ceramics Analysis and Reliability Evaluation (SCARE) design program. 


INTRODUCTION 

The fracture behavior of ceramics differs considerably from that of metals 
This is due to several factors: Little or no plastic deformation exists at 

crack tips in ceramics. They also possess lower fracture toughness, and sub- 
critical cracks grow under static or dynamic stresses, even in ambient environ- 
ments. Ceramics are also not generally homogeneous. They may exhibit one or 
more types of defects with different distribution of sizes. The observed varia 
bility in ceramic strength is due to the presence of microcracks, the large 
range of flaw sizes, and the residual stresses arising from anisotropic contrac 
tion upon cooling. The design of load-bearing components using structural cer- 
amics requires a complete understanding of the statistical nature of fracture 
in brittle materi al s . 

The reliability and consistency of design strength data for ceramics is 
relatively poor when compared to the design data for materials like metals and 
plastics. The reason for this is that ceramics contain a large number of 
microcracks with various sizes and orientations which, along with residual 



stresses from processing, cause a wide scatter in strength measurements. The 
great variation in strength of nominally identical specimens necessitates test- 
ing large numbers of samples to fully describe the failure phenomena. However, 
this poses a dilemma since the costs associated with testing equipment, speci- 
men preparation, and precautions to ensure against damage of specimens are high. 
Thus, it is desirable to test an optimum, but relatively small, number of sam- 
ples and to develop analytical procedures that yield material parameters with a 
high degree of accuracy. 

A number of commonly used uniaxial 1 y loaded test specimens can be used to 
measure ceramic fast-fracture strength over a wide range of temperatures. 

These material characterization tests are derived from observed failure modes. 
The Army Specification (ref. 1) describes a test method for determining the 
flexure strength, or modulus of rupture (MOR), of brittle beam test specimens 
at ambient temperatures. Because of the scatter in observed strength values, 
many samples (typically 30 or more MOR bars for a specified temperature) are 
required. Fractography is required after fracture to separate flaw populations 
and to identify failure modes. The major advantages of this testing method 
include simple fixtures, low cost, and the ability to cut small test specimens 
out of component bulk. 

However, an inherent weakness of flexure testing is that it emphasizes sur- 
face flaws relative to volume flaws. It also needs a greater order of magni- 
tude of volume and area scaling between the specimen and the entire component 
geometry. To alleviate these problems, considerable effort has been focused on 
developing and using tensile fast-fracture tests. One form of these tests, 
recently developed at Oak Ridge National Laboratory (ORNL), uses a hydraulic 
self-aligning grip system (ref. 2). With the tensile test specimen, the entire 
gage material volume (which is generally much larger than in the flexure speci- 
men) is loaded to the same maximum stress, and many successful tests produced 
near zero bending during monotonic tensile and cyclic tension-tension loading. 
With the emerging importance of ceramic tensile testing, the need for a smaller 
number of specimens, more advanced statistics for data reduction, and more 
selective testing than with MOR bars is evident. 

Two main classes of fracture theory have been previously discussed in the 
literature. The first is derived from Griffith's flaw theory (refs. 3 and 4) 
which assumes the presence of cracks of specified sizes, shapes, and arbitrary 
orientations. There is always one crack with the least favorable orientation 
that causes failure. The second class consists of the weakest link theory (WLT) 
such as that of Weibull (ref. 5). This theory is purely statistical in nature. 
Failure in ceramics is due to a single weakest flaw where the local stress 
reaches a critical value. The inability of current inspection procedures to 
fully characterize flaws and the observed size effect on strength makes it nec- 
essary to use statistics and reliability analysis to predict failure of brittle 
material components. 

Weibull (ref. 6) used the WLT to develop a probabilistic approach to _ _ 

account for the scatter in the fracture strength and the size effect of brittle 
materials. He assumed a unique cumulative strength distribution for the uniax- 
ial fracture data obtained from simple test specimens. Weibull also adopted 
uniaxial fracture test parameters for calculating the failure probability in 
multidimensional stress states. The concept focuses on calculating the risk of 
rupture by averaging the normal tensile stress at any point in principal stress 
space. This approach is plausible, but rather arbitrary, and requires numeri- 
cal modeling. Consequently, other models were proposed, and the most widely 
used model assumes that the principal stresses act Independently (PIA) 
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(ref. 7). However, the PIA hypothesis may lead to unsafe estimates of failure 
probability because it neglects the effects of the interaction of combined 
local principal stresses (ref. 8), 

Ceramic strength is well represented by the Weibull distribution which is 
fully characterized by its statistical parameters. Batdorf and Crose (ref. 9) 
proposed a new physically based statistical theory, which describes the mate- 
rial strength by a flaw-density function. This strength parameter is usually 
expressed in terms of a power function with appropriate flaw-density constants. 
It is then a simple matter to derive relationships between the Weibull parame- 
ters and the Batdorf flaw-density constants, so that the same material strength 
distributions are consistently obtained from both. Rufin et al . (ref. 10) 
further extended the Batdorf statistical theories to investigate effects of dif- 
ferent flaw geometries and mixed-mode fracture criteria on overall material 
reliability. Gyekenyesi and Nemeth (refs. 11 and 12) developed a public domain 
computer program incorporating the preceding theories which could be coupled to 
a general purpose finite element code such as MSC/NASTRAN (ref. 13) to predict 
the fast-fracture failure probability of ceramic components due to the presence 
of volume- and surface-type flaws. 

The primary objective of this paper is to describe the estimation of the 
Weibull strength parameters and the Batdorf flaw-density constants for volume- 
and surface-flaw-induced fracture in ceramics. An extensive literature survey 
of various statistical estimation techniques was conducted and discussions were 
held with numerous researchers in the field. The goal was to quantify the 
uncertainty by using analytical methods. Many modern Weibull parameter estima- 
tors (ref. 14) are efficient and provide confidence intervals. Mann et al . 

(ref. 15) discuss moment estimators, least-square estimators, best least invari- 
ant estimators, and maximum likelihood estimators (MLE's). The evaluation of 
shape and scale parameters of the two-parameter Weibull distribution is 
explained initially in this paper by using the least-squares analysis (ref. 16). 
This method obtains, analytically, a unique line of best fit to describe the 
scatter in ceramic strength. 

The MLE's have the smallest mean-square error, and Bain and Antle 
(ref. 17) found them to be the most efficient of the available methods. How- 
ever, MLE's are highly biased for small samples, and, consequently, the MLE of 
the shape parameter is often unbiased to eliminate the deviation between the 
sample and true population. Important design criteria should not be derived 
from small sample sizes because uncertainty is known to increase as the number 
of samples decreases. Menon (ref. 18) describes the derivation of the maximum 
likelihood equations based on the assumption that the location parameter 
(threshold stress) is known. This method was further extended by Cohen 
(ref. 19) for singly censored samples. For small to moderately large sample 
sizes, distribution percentiles of the MLE's and unbiasing factors for the 
shape parameters have been computed and tabulated by Thoman et al . (ref. 20) 
through the use of Monte Carlo simulation procedures. Jeryan (ref. 21) devel- 
oped a FORTRAN computer program to calculate the MLE's and other related 
statistical parameters for uncensored statistics. This paper summarizes the 
equations and procedures that are necessary for evaluating the MLE's of the 
shape and scale parameters, the Weibull mean values, and the corresponding 
standard deviations (ref. 22) by using theories of concurrent flaw distribu- 
tions (ref. 23). The Weibull log likelihood equations for censored statistics 
were developed by Nelson (ref. 24) and used to calculate the MLE's for the cen- 
sored data using Newton-Raphson iterative techniques. Procedures to evaluate 
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90-percent confidence intervals for MLE's of the shape and scale parameters and 
the unbiased estimates of the shape parameters for uncensored statistics are 
explained. No rigorous studies are known to the authors that evaluate the con- 
fidence intervals for MLE's of shape and scale parameters and the unbiased 
estimate of shape parameter for censored statistics. However, a Monte Carlo 
simulation is developed and tested for censored statistics, and the results are 
included in this paper. 

An outlier is a data value that is far from the rest of the sample. There- 
fore, techniques are needed that detect outliers to reveal faulty products as 
well as to improve poor test data. Stefansky (ref. 25) describes a procedure 
for calculating critical values of the maximum normed residuals while testing 
for outliers. This technique is based upon the normal distribution, and its 
application to Weibull data is of concern. However, it was used by Neal et al . 
(ref. 26) to calculate the critical values for different significance levels. 
Sample size-dependent polynomial approximating functions have been developed to 
compute the critical values at 1-, 5-, and 10-percent significance levels, and 
these functions are described in this paper. 

Many goodness-of-f i t tests based on empirical distribution functions (EDF) 
have been extensively discussed by D'Agostino and Stephens (ref. 27). The EDF 
statistics measure the discrepancy between the EDF and a given distribution 
function. Five of the leading EDF statistics, including the Kol mogorov-Smi rnov , 
the Cramer-von Mises, the Kuiper, the Watson, and the Anderson-Darling have 
been examined by Stephens (ref. 28). The Kol mogorov-Smi rnov (K-S) test for 
goodness-of-f i t is the most widely used in comparing an empirical distribution 
function with the distribution function of the hypothesized distribution for 
any finite sample sizes (ref. 29). Critical values of the maximum absolute dif- 
ference between sample and population cumulative distributions have been tabu- 
lated by Massey (ref. 30). Percentage points in Kol mogorov-Smi rnov statistics 
for probabilities from 0.01 through 0.20 were tabulated by Miller (ref. 31) for 
different sample sizes. This work was further extended by Amstadter (ref. 32) 
in evaluating values for probabilities from 0.30 through 0.99. The Anderson- 
Darling (A-D) statistic A 2 (ref. 33) has proven more sensitive to the discrep- 
ancies in the tail regions between the EDF and the cumulative distribution 
function. Lewis (ref. 34) discusses the Anderson-Darling statistic for a fully 
specified distribution. This paper summarizes the steps to calculate the 
Kol mogorov-Smi rnov and the Anderson-Darling goodness-of-f i t statistics by using 
polynomial approximating functions to represent available tabular values from 
the literature. 

Confidence bands (ref. 35) on the Weibull line can be calculated and plot- 
ted (ref. 36) for different confidence levels by making use of the mean and the 
standard deviation of the sample. This paper also explains the calculation of 
90-percent confidence bands about the Weibull distribution. 

Also, several currently used theories that can be used for calculating the 
Batdorf flaw-density constants from the Weibull distribution statistical param- 
eters are selectively included in this paper. The most widely used model among 
these theories assumes that the cracks are shear-insensitive. Different rela- 
tionships are obtained, however, if shear-sensitivity of the flaws is included. 
Finally, all the statistical techniques summarized here were verified with sev- 
eral example problems from the open literature, including the problem of deep- 
groove ball bearings (ref. 37) and the recently obtained Elecktroschmel zwerk 
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Kempten (ESK) hot isostatically pressed HIPped silicon carbide (ref. 38) 
strength data. 


PROGRAM CAPABILITY AND DESCRIPTION 

> 

Volume- and surface-flaw-based reliability analysis has been implemented 
in both the SCARE1 and SCARE2 versions of the postprocessor program (refs. 11 
and 12). The SCARE 1 version of the code uses only elemental centroidal princi- 
pal stresses to calculate the reliability for both volume and surface elements. 
In SCARE2 , all linear or quadratic QUAD8 shell elements are further discretized 
into nine subelements with interpolated centroidal principal stresses and sub- 
element areas used to perform reliability analysis. A similar discretization 
of brick volume elements into 27 subelements is performed for volume-flaw- 
based reliability studies. Both versions of SCARE have broad capabilities: 
Temperature-dependent multiple-flaw populations, statistical material parame- 
ters, several crack configurations, and a number of well-known fracture criteria 
can be specified. The existing architecture and basic computational elements 
of the programs were described in previous publications (refs. 11 and 12). Fig- 
ure 1 shows the flowchart for calculation of shape and scale parameters of the 
two-parameter Weibull distribution using the least-squares analysis and the max- 
imum likelihood methods for volume- and surface-flaw-induced fracture in ceram- 
ics with complete and censored samples, and the calculation of many other 
related statistical quantities. For convenience, these latest additions to the 
SCARE program can also be used independently of the MSC/NASTRAN postprocessing 
or without performing a reliability analysis. 


INPUT INFORMATION 

SCARE input requirements are grouped into three major categories: the 

Master Control Deck, the Specimen Deck, and the Structures Deck. The details 
of the deck functions are explained in reference 11. For the calculation of 
shape and scale parameters of the two-parameter Weibull distribution, the data, 
which includes the fracture origins of the specimen, is entered in the Specimen 
Deck. Initially, for the statistical analysis of both complete and censored 
samples, each specimen fracture mode is observed by fractography or other means 
and identified as failed either by volume-flaw- or surface-flaw-induced frac- 
ture, V and S, respectively. If there is any uncertainty regarding a failure 
mode, it is input as an unknown failure mode U. The user can select any of 
the following three options of calculations for complete samples: Volume flaw 

analysis only, surface flaw analysis only, or both in one execution. For the 
last selection, two separate extreme fiber fracture stress data sets (which may 
be completely identical) are required for the analysis. The same user options 
are available for censored samples. When the third choice is selected, volume- 
flaw analysis will be done by treating surface-flaw and unknown-flaw failure 
mode data as suspended items; surface-flaw analysis will be done by treating 
the volume-flaw and unknown-mode data as censored information. 


OUTPUT INFORMATION 

The details of SCARE output for a complete MSC/NASTRAN finite element 
analysis are discussed in previous publications (refs. 11 and 12). The output 
pertaining to the calculation of the Weibull strength parameters identifies the 
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method of solution, the control index used for experimental data, the number of 
specimens in each batch, and the temperature of each test. In addition, the 
output echoes the input values of all extreme fiber fracture stresses with 
proper failure mode identification. Any data value that is far from the rest 
of the sample is detected as an outlier, and its corresponding significance 
level (1, 5, and 10 percent) is printed. Furthermore, the biased and the 
unbiased values of the shape parameter, the specimen characteristic strength, 
the upper and lower bound values at 90-percent confidence level for both the 
shape parameter and the specimen characteristic strength, the specimen Weibull 
mean value, and the corresponding standard deviation are printed for each speci- 
fied temperature. For censored statistics, the above values are generated 
first for volume-flaw analysis and, subsequently, for the surface-flaw analysis. 

The Kolmogorov-Smirnov goodness-of-f i t test is done for each data point, 
and the corresponding K-S statistic factors D + and D - are listed. Similarly, 
the K-S statistic for the overall population is printed along with the percent- 
age significance level. This overall statistic is the absolute maximum of indi- 
vidual data D + and D~ factors. For the Anderson-Darling goodness-of-f i t test, 
the A-D statistic A 2 is determined for the overall population, and its associ- 
ated significance level is printed. 

The next table of the output contains fracture stress data, the correspond- 
ing Weibull probability of failure values, the 90-percent upper- and lower- 
confidence band values about the Weibull line, and the median rank value for 
each data point. These statistical quantities are calculated by using either 
tabular values or approximating polynomial functions which are now incorporated 
into the SCARE program. Finally, the material parameters used In component 
reliability calculations are listed as a function of temperature. These 
include the biased Weibull modulus, the Batdorf crack-density coefficient, and 
the material Weibull scale parameter or unit volume and area characteristic 
strength, whichever is appropriate. 


THEORY 

The weakest link theory of material strength considers not only the size 
effect and loading system, but also the variation in fracture strength due to 
intrinsic defect distributions. Paluszny and Wu (ref. 39) discuss evidence of 
adequate correlation between the reliability predictions of complex structures 
based on simple test bars using the Weibull methods and the experimental 
results. It is possible to use the Weibull distribution to model a broad range 
of instantaneous failure rates and also to represent several component failure 
mechanisms. It has been shown that experimental fracture strength data 
obtained from uniaxial 1 y loaded specimens is best approximated by the three- 
parameter Weibull distribution, defined by 


P 


f “ 



( 1 ) 


where Pf Is the probability of failure, a is the applied tensile stress, a u 
is the threshold stress (location parameter usually taken as zero for ceramics), 
ct 0 is the scale parameter with dimensions of stress x (volume) 1/m , m is the 
shape parameter (or the Weibull modulus) which measures the degree of strength 
variability, and V is the stressed volume. A similar equation can be written 
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for area-flaw-induced failure with area parameters replacing corresponding 
volume variables. 

The location parameter is the value of a below which the failure proba- 
bility is zero. Similarly, the scale parameter corresponds to a stress level 
where 63.2 percent of specimens with unit volumes would fracture. The Weibull 
modulus is a measure of the strength dispersion, and as m increases, the 
strength becomes increasingly deterministic. Values of m = 3.44 and m = 1 
give, respectively, the normal (approximate) and true exponential distribution 
curves as shown in figure 2. For large values of m (m > 40), such as those 
obtained for ductile metals, the scale parameter corresponds to the material 
ultimate strength. Typically, for brittle materials, the Weibull parameters 
are determined from simple geometry and loading conditions, such as beams under 
flexure or cylindrical or flat specimens under uniform uniaxial tension. 


The Outl i ers 


Before computing the estimates of the Weibull parameters, the available 
specimen fracture stress data must be carefully examined for outliers. Very 
often a data set may contain one or more values at the extremes which may not 
belong to the main trend of the overall population. Such data points or values 
are labeled as outliers. It is not always possible to notice these outliers 
without having a complete and thorough knowledge of the manner in which the 
data is obtained. However, from statistical computations, the outliers can be 
detected at different significance levels. 

At the start of the outlier calculations, the sample mean of and the sam- 
ple standard deviation s are calculated. From these values, the normed resid- 
uals q for each sample are obtained as follows: 


r 


i 



i = 1,2, . 


( 2 ) 


where N is the sample size. The normed residuals are normalized deviations 
of the data from the "center" of the data. The absolute maximum of the normed 
residual (MNR statistic) is compared to the critical value (CV), associated 
with the sample size. Critical values are calculated from the following 
equation: 


CV = \ ^ j (3) 

y/H T N - 2 + t 

where t is the (1 - a/2N) quantile of the t-di stri bution with N - 2 degrees 
of freedom, and a is the significance level. Separate polynomial approximat- 
ing functions for sample sizes to 30 and above 30 have been developed for calcu- 
lating CV at the appropriate significance level and are included in the SCARE 
program. These functions evaluate the critical values at 1-, 5-, and 10-percent 
significance levels. If the MNR statistic is smaller than the three critical 
values, then no outliers are detected. However, if the MNR is larger than at 
least one of the three critical values, the corresponding data value with the 
MNR statistic is detected as an outlier with the appropriate significance level. 
Once all the deviant points are detected, each potential outlier is retested 
against the remaining "good" data, and the results of this test are printed in 
the program output. 
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It Is sometimes appropriate to omit the outliers and continue with the 
remaining analysis treating the reduced sample as a new sample or as a censored 
sample. But the rejection of outliers, strictly based on statistical analysis, 
can be misleading. Therefore, this technique of detecting the outliers should 
be used as a means of ascertaining the relative validity of the available data. 
Figure 3 shows the existence of a potential outlier in a data set for sintered 
Si 3N4 generated at NASA Lewis Research Center in cooperation with Oak Ridge 
National Laboratory. It appears that competing failure modes exist, and an 
outlier point is seen at the high-strength range of the data. Since fractogra- 
phy was not performed, the two Weibull distributions shown are assumed. 

Weibull parameters were only calculated for the combined distribution, which 
was treated as a complete sample, and its plot is also shown. 


Estimation of Statistical Material Parameters 


The flexural test failure probability can be expressed in terms of the 
extreme fiber fracture stress of, or MOR, using the Weibull form as 


P 


f ~ 


_ m 
exp -Cof 


- exp -C(M 0 R) m ] 


( 4 ) 


where C = (I/oq) 111 and 09 is the volume/area specimen characteristic 
strength, or characteristic modulus of rupture MORp. In equation ( 1 ), o 0 is 
based on unit volume; whereas in equation ( 4 ), 09 includes effects of the spec- 
imen volume. A sample may only contain data from specimens that failed either 
from a surface flaw or from a volume flaw. Such a sample is usually called a 
complete sample. In some cases, specimens in a sample may not fail by one fail- 
ure mode only, but by several competing failure modes. Such a situation exists 
in the case of concurrent flaw populations. Each specimen fracture origin is 
identified by fractography and then labelled as failed either due to a volume 
flaw or a surface flaw. For the computer input, these failure modes are denoted 
as V and S, respectively. Occasionally, some specimen fractures may not be 
identified directly due to uncertainty in the dominance of failure modes. Such 
data are grouped into an uncertain category and are denoted by U. Thus, the 
sample containing data from the mixture of failure modes, and sometimes inclu- 
sive of unknown failure modes, is referred to as a censored sample. 


Many methdds~are available to estimate the statistical material parameters 
from experimental data. The success of the statistical approach depends on how 
well the probability function fits the actual material measurements, especially 
how accurate the value of m is. The two most popular models used to evaluate 
C and m are the least-squares analysis and the maximum likelihood methods. 
These statistical techniques are described in detail in the succeeding sections. 


The Least-Squares Analysis 

The survival probability of a ceramic specimen in flexure strength testing 
is defined as 

P $ = exp[- C of] (5) 
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This equation can be linearized by taking the natural logarithm twice to obtain 


Sin 



an 


an 



an C + m an 


( 6 ) 


If the failure probability Pf is determined from conducting N tests and num- 
bering the observed fracture stresses as 01 , 02 , ■ • oj in an 

ascending order, then 


W - NTT 


(7) 


where i is the median rank. A more sophisticated treatment using median rank 
regression analysis gives the following result (ref. 36): 


Pf (°i > = 


i - 0.3 
N + 0.4 


( 8 ) 


The difference between equations (7) and (8) becomes insignificant for large 
values of N. However, equation (8) is used in the remaining analysis as well 
as in the code. For the least-squares analysis, it is necessary to obtain the 
line of best fit with slope m and an intercept b which, as can be noted 
from equation (6), is equal to the natural log of C. By differentiating the 
sum of the squared residuals with respect to m and b separately, setting 
the results equal to zero, and then simultaneously solving for m and b, we 
obtain 


- gViE«i 


"2>? - X>(2>i 

where y^ and are also obtained from equation (6). Thus, the material 
constants are calculated from known fracture stress data afj at a specified 
temperature. With censored data, one cannot directly use the median rank 
regression analysis as given in equation (8) because of the presence of compet- 
ing failure modes. To take into account the potential influence of the sus- 
pended items, Johnson (ref. 40) developed the rank increment technique, which 
is calculated from the following equation taken from reference 36: 


Rank increment = 


(N + 1) - (previous adjusted rank) 

1 + (number of items beyond present suspended item) 


( 11 ) 


As mentioned earlier, all observed fracture stresses are arranged in 
ascending order, and rank increment values are calculated for each failure data 
by using equation (11). For volume-flaw analysis, all MOR's designated as 
V's are considered as failure data; for surface-flaw analysis, the S's are 
considered as failure data. According to the procedure discussed in refer- 
ence 36, the new adjusted rank values are obtained by adding the rank increment 
value to the previously adjusted rank. These adjusted rank values are then 
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used to calculate the failure probability Pf by using median rank regression 
analysis (i .e. , eq. (8)) . Finally, the Weibull parameters m and C are 
obtained by solving equations (9) and (10), except that N now only includes 
the number of failed samples and not the total samples as in the case of com- 
plete (uncensored) data. 

The least-squares analysis is a particular case of the maximum likelihood 
method, if we assume that the error is normally distributed with zero mean and 
constant variance. Somerville and Bean (ref. 41) compared the least-squares 
analysis and the maximum likelihood methods for different sets of conditions. 

The study concluded that when the data fits the Weibull distribution, both tech- 
niques produce substantially the same results. The least-squares analysis is 
most useful in situations where the underlying probability distribution is not 
clearly established, and it is more robust for minor departures from the 
Weibull distribution. In addition, the least-squares analysis may give better 
.results for censored samples. While discussing the statistical approach to 
engineering design in ceramics, Davies (ref. 42) noted that the least-squares 
analysis gives better estimates of the parameter if a single anomalous high- 
strength value is present in the sample. King (ref. 43) determined the accuracy 
of the estimates of the Weibull parameters by using the prediction interval 
method, and he showed that the sampling distribution of the regression coeffici- 
ents follows a t-di stri bution . This is true only if the error is normally dis- 
tributed with zero mean and constant variance. Therefore, the least-squares 
analysis method is usually not suited to calculate confidence intervals and 
unbiasing factors and to take into account any uncertainty in the available 
data. For that reason, the maximum likelihood method has been extensively used 
in Weibull analysis and is discussed in detail here. 


The Maximum Likelihood Method 

The maximum likelihood method is often used because of certain inherent 
properties. The likelihood equation from which the MLE's are obtained has a 
unique solution. In addition, as the sample size increases, the solution con- 
verges to the true values of the parameters, and, therefore, the solution is an 
asymptotically normal and asymptotically efficient estimate of the parameters. 
Jeryan (ref. 21) compared the maximum likelihood method to the method of 
moments and the pseudo-least-squares analysis for random sample sizes with 
assigned values for the Weibull parameters. It was concluded that the MLE of 
the Weibull slope showed the least bias and variance. According to Sonderman 
et al. (ref. 44), the maximum likelihood method is a wholly parametric tech- 
nique, dependent on the form of the Weibull distribution to obtain its esti- 
mates. A comparative study (ref. 44) between the mean order, the hazard 
plotting, and the maximum likelihood estimates indicated that the MLE of the 
shape parameter was reasonably accurate for small samples. For complete as 
well as censored samples, no ranking functions or linear regression analysis 
are used while using the maximum likelihood method. 

The likelihood equation is written by expressing the probability of obtain- 
ing the observed data. For the complete sample, the likelihood function is 
defined by 


N 

L = II f(o f .) = f(a f] )f(o f2 ) . . . f(o f|y| ) (12) 

1 = 1 H 
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In the analysis of failure of brittle materials subjected to uniaxial stress 
states, the Weibull probability density function is 


f (o f > 



and the corresponding likelihood function is 



(13) 


(14) 


In equation (14), the "likelihood" of the sample failure data af ] , a f 2 , . . ., 
CTf i , . . afN is a function of the Weibull parameters m and oq. The val- 
ues of m and oq which maximize the likelihood function are determined by 
the maximum likelihood method. By differentiating the logarithm of the likeli- 
hood function with respect to m and oq, and then equating the resulting 
expressions to zero, we obtain the MLE's of m and oq by solving the result- 
ing simultaneous equations. 


In the complete sample case, the MLE of the shape parameter m, denoted by 
m, should satisfy the following equation (ref. 36): 


/ 


N 

y~l (a f . ) m Q,n(a f . ) 
i=l 71 71 



N 


i = l 



(15) 


The Newton-Raphson iterative technique is used to solve equation (15) and to 
calculate m. Furthermore, the MLE of the Weibull parameter oq, denoted by 
oq, is evaluated by solving the following equation: 



However, if the sample is censored and the effect of a competing failure mode 
is taken into account, the following equation is solved by iterative techniques 
to obtain the value of m (ref. 36): 


N 

£ (o f .An(o f .) 

i=l 71 71 

N 
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I 
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i = l 


&n(cr f .) - l = 0 
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(17) 
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In the preceding equation, r is the remaining number of specimens failed 
either by volume flaw or by surface flaw only. Finally, ctq is obtained from 


= 


e 



( 18 ) 


In the previous cases, the MLE of the Weibull parameter C, as used in equa- 
tion (4), is calculated by using the previously explained conversion technique. 
Incidental ly.^the MLE's of m and 09 for the volume-flaw analysis are denoted 
by my and oq\j\ for the surface-flaw analysis by me; and 093, respectively. 
Similarly, the parameter C is denoted as Cy and C$. 

As shown in figure 1 , the least-squares analysis is always done first, 
irrespective of the method of analysis chosen. However, if the maximum likeli- 
hood method is selected, the initial estimate of the shape parameter for the 
iterative technique is set equal to the value of the parameter that is obtained 
from the least-squares analysis. This is convenient for computation, and the 
number of iterations in the Newton-Raphson technique is arbitrarily limited to 
50 . If no solution is found after 50 iterations, the maximum likelihood method 
is terminated, and the results from the least-squares analysis are used for fur- 
ther calculations in the program for both complete and censored samples. 


Unbiasing Factors and Confidence Intervals 

In the previously described maximum likelihood method, the MLE of the 
shape parameter is a biased estimate that depends on the number of specimens in 
the sample. Consequently, this parameter may not represent the true popula- 
tion. To minimize the deviation between the sample and the total population, 
the shape parameter estimate should be unbiased. Statistically, an estimator 
of any parameter is called an unbiased estimator if the mean of the estimators 
of its sampling distribution equals the value of the parameter itself. The non- 
zero difference between the estimator of the parameter and the parameter Itself 
is called the bias of the estimator. Estimates are random variables, and a 
thorough understanding of the accuracy of the parameters estimation is needed. 
This is usually achieved by obtaining the confidence intervals of the distribu- 
tion parameters. If the variable 6 and a quantity 0*. respectively, define 
the estimated and the true value of the parameter, then according to Siddall 
(ref. 45), the confidence interval is defined as an interval in 0 that will 
include the unknown value 0* with a probability [J, where |3 is called the 
confidence level or confidence coefficient. The related uncertainty with the 
parameter estimate is based on a sample of data, and a different data set of 
the same size will yield different results. In general, the confidence inter- 
vals depend upon the sample data and not upon the parameters themselves. Conse- 
quently, the unbiasing factors, as well as the confidence intervals, play a 
major role in the analysis. 

Thoman et al . (ref. 20 ) obtained the unbiasing factors for the MLE of m 
for complete samples by initially obtaining the MLE of m for a standard expo- 
nential density function, which is a two-parameter Weibull density function as 
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given in equation (13) with unit values of m and uq . The unbiasing factors 
for the MLE of m in reference 20 depend on sample size. The unbiased esti- 
mate of m is obtained by multiplying the biased estimate of m with the 
unbiasing factors. The percentage points of the distributions of MLE's were 
evaluated and then used to calculate the confidence intervals for MLE's at sev- 
eral confidence levels for complete samples (ref. 20). McLean and Fisher 
(ref. 46) discuss the uncertain nature of Weibull parameters estimations. 
Figures 4 and 5, taken from reference 1, show 90-percent confidence bounds for 
the Weibull modulus and the characteristic strength. As can be seen from the 
graphs, the error or the uncertainty in estimates from a small sample size is 
very large. The curves in figures 4 and 5 are obtained from tables 1 and 4 of 
reference 20. 

A method of obtaining definite confidence intervals for censored samples 
has not yet been developed because of the additional complexity of competing 
failure modes. However, the effects of competing failure modes were examined 
through a Monte Carlo simulation on a limited trial basis by the authors. The 
investigation was performed by assigning two specified Weibull distributions as 
the parent populations, running a Monte Carlo simulation to mimic the competing 
failure modes, and recovering the MLE's for those trials, where a predetermined 
number of failures from a given mode occurred. It is estimated that 10 000 
trials provide reasonable accuracy when calculating the 90-percent confidence 
limits of the parameters. Wetherhold (ref. 47) discusses in detail the statis- 
tical reasoning for the proper selection of the number of trials in the Monte 
Carlo simulation. With this method, the effect of sample size, including a 
given percentage of failure modes from the trial runs, indicated that the per- 
centage points of the distribution asymptotically approached the complete sam- 
ple case and depended on the number of failures for each failure mode. 

It is clear from the preceding discussion that a knowledge of the upper 
and lower bound values of the MLE's of both parameters at any confidence level 
is necessary. This enables the user to estimate the uncertainty in the parame- 
ters as a function of the number of specimens. It is general practice to 
obtain these bounds at 90-percent confidence level, and, therefore, 5 and 
95 percentage points of distribution of the MLE's of the parameters have been 
incorporated into the SCARE program, with data taken from reference 20. Fur- 
thermore, this also assures that if different data sets are used in the maximum 
likelihood method, then 90 percent of the time, the MLE's of the parameters 
will lie between the upper and lower bound values. The following equations 
describe the manner in which the different limit values are calculated at 
90-percent confidence level. 
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where a minus sign has been used in equation (21) to account for the data for- 
mat. To calculate the upper and lower bound values, 5 and 95 percentage points 
of the distributions of the MLE's are used, respectively. All these values 
depend on the biased value of the MLE of the shape parameter, denoted by 
mBiAS ■ Jeryan (ref. 21) used the unbiased value of the shape parameter MLE to 
calculate the upper and lower bounds, while Thoman et al . (ref. 20) used the 
biased value as done here. Different approximating equations for these bounds 
can be found in reference 36 for complete samples. 


The Weibull Mean and the Standard Deviation 


After evaluating the MLE's of the Weibull parameters, the Weibull mean y 
and the standard deviation s can be obtained by using the following two 
equations (ref. 22): 


V- = 




/ 






1/2 


(23) 


(24) 


where T denotes the gamma function. These equations are valid for both area- 
and volume-based flaw strength. Equations (23) and (24) indicate that as m 
increases, the mean of the distribution approaches the specimen characteristic 
strength, and the standard deviation approaches zero. We can also note that 
the specimen characteristic strength ctq is related to the material scale 
parameter a 0 through the specimen effective volume/area relationship of 

°ov = ^ev v e /mv or °os = °es A e /nis - 


Goodness-of-Fi t Tests 

Graphical techniques can often be used to determine whether or not a set 
of observations belongs to a population distribution. A discussion of these 
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methods for the Weibull and the lognormal distributions is given in refer- 
ence 15. Subjective judgment is needed to test the goodness-of-f i t of the data 
to the assumed distribution. In many instances, it is difficult to decide if 
the hypothesized distribution is valid. Therefore, classical statistics tests 
are used in verifying the available data. 

In general, a statistic is a numerical value determined from a sample by 
using statistical computations. The difference between an empirical distribu- 
tion function (EDF) and a given distribution function is called an EDF statis- 
tic. There are two major classes of statistics, and they differ in the manner 
in which the functional (vertical) difference between the two distribution func- 
tions is calculated. The supremum class statistic calculates the largest and 
smallest vertical differences between the two distribution functions. On the 
other hand, the quadratic class statistic measures the discrepancy between the 
two distributions through squared differences and the use of an appropriate 
weighting function. 

The Kolmogorov-Smi rnov goodness-of-f i t statistic belongs to the supremum 
class and is very effective for small samples. To begin the K-S test, the sam- 
ple is arranged in ascending order, and the EOF, F^fx), is obtained by using 
the following expressions (ref. 27): 


F N <x> = 0, 

x < X <D 




F n ( x) - 

X ( i ) - x ( X (x»l>' 

i = 1,2, . . 


> (25) 

l 

F n (x) = 1, 

X ( N) i * 


1 

i 
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where x] , X 2 x^ are from the sample. Figure 6 shows the step func- 

tion calculated from the data. It Is obvious that as x increases, the step 
function F(yj(x) steps up as an individual sample observation is reached. For 
any x, F^(x) denotes the proportion of observations less than or equal to x; 
F(x), defined as the distribution function, is the probability of an observa- 
tion less than or equal to x. The main intent of the K-S test is to calculate 
the well-known EDF statistic D, which is the maximum absolute vertical differ- 


ence between the hypothesized cumulative distribution function F(x) and the 
EDF function F^( x ) . As shown in figure 6, the statistic D is obtained by 
initially evaluating two other statistics D + and D - , the largest vertical 
differences when F^(x) is greater than F(x) and the largest vertical differ- 
ences when Ffj(x) is smaller than F(x), respectively. In the SCARE program, 
all three statistics are calculated by using the following expressions 


(ref. 28): 




(26) 


For ceramics design, the F<x)}'s are equal to Pf's and are calculated 
by using equation (4). These are also referred to as the predicted failure 
probability for each fracture stress level and are evaluated by using the 
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biased values of the MLE's of the Weibull parameters. In the censored samples, 
the presence of two concurrent flaw populations makes it essential to calculate 
the cumulative survival probability first. This is accomplished by using 
equation (5) separately for volume and surface flaws, and then by multiplying 
the survival probabilities of the two concurrent flaw populations. Finally, the 
resultant predicted failure probability at each stress level is calculated by 
using the entire sample. 

It is customary to express the EDF statistic D with an associated signif- 
icance level a, which is the probability of obtaining a higher maximum absolute 
difference assuming that the data follows the theoretical distribution. This 
probability is a function of D and the sample size N. Miller (ref. 31) 
described a general formulation for evaluating a as follows: 


a = 


exp 


2D 


('-¥)] 


(27) 


For sample sizes of 30 or less, the values of a obtained by using equation 
(27) are equal to the values given in reference 32. However, for sample sizes 
greater than 30, a separate formulation is used based on the discussion by 
Brunk (ref. 48), and is given by I 



In the SCARE program, the EDF statistics D + or D~ are calculated for every 
specimen fracture stress level. Then the EDF statistic D and the associated 
significance level are determined for the total sample size. A higher value 
of a indicates that the data fits the proposed distribution to a greater 
extent. 

The Anderson-Darling statistic A 2 belongs to the quadratic class and is 
a powerful statistic for goodness-of-f i t measure. This approach measures the 
discrepancy between the EDF and the hypothesized cumulative distribution func- 
tion by using the Cramer-von Mises formulation (ref. 27): 


Q = 



{F N ( x) - F(x>} 2 i|/(x)dF(x) 


(29) 


where y(x) is a suitable weighting function applied to the squared difference. 
Anderson and Darling (ref. 33) introduced the A 2 statistic by defining y(x) 

= [ {F< x ) } { 1 - F(x)}]- 1 in equation (29). This statistic is more suitable for 
identifying the discrepancies in the tail regions of the distribution. The con- 
cept of the probability integral transformation (PIT) is discussed in detail in 
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reference 27. In this procedure, a new random variable Z, which is uniformly 
distributed between 0 and 1, is introduced so that Z = F(X), when F(x) is the 
true distribution of X. The Z^'s are then arranged in ascending order and 
used to calculate the statistic which is given by (ref. 27): 


A 


2 



1) an{Z (1) } + dn{l 


Z (N+l-i 


(30) 


In this case, Z^'s are the predicted failure probabilities obtained from equa- 
tion (4). It can be shown that for N > 3, the distribution of A^ is suffici 
ently independent of the sample size. 

Lewis (ref. 34) discussed the asymptotic distribution and gave the values 
for the A^ statistic and the corresponding significance levels. These values 
are summarized in reference 27. We have developed several approximating polyno 
mial functions to interpolate these values. If the A^ statistic is too large, 
the data should be rejected because the EDF fits the proposed distribution 
poorly. Otherwise, the significance level for the sample is calculated by 
using the A^ statistic and the corresponding polynomial function. Once 
again, the higher the value of a, the better the data fits the proposed 
di str i bution . 


Confidence Bands on the Weibull Line 

Kanofsky and Srinivasan (ref. 35) described a method for constructing a 
confidence band on a continuous distribution function which may have normal, 
exponential, and many other functional forms. They introduced and defined the 
random variable L, (the statistic L) as 


L = 


max 

x 



VV 


- F(x; x,s) 


(31) 


where p p , E p , x, and s are, respectively, the population mean, the popula- 
tion standard deviation, the sample mean, and the sample standard deviation. 

The population parameters p p and E p are unknown values. The sample is 
drawn_from a continuous population distribution. In equation (31), the statis- 
tic L is the maximum absolute difference between the two cumulative distribu- 
tion functions. The distribution of the statistic L is obtained in two 
stages: (1) The closed-form expressions for L are derived. These expres- 

sions are functions of x and s; and (2) these expressions derive the distri- 
bution of L by using Monte Carlo simulation. The random variable L is then 
used to construct different percentage confidence bands on F ( x ; p p , E p ) cen- 
tered around F(x; x, s). This is similar to the use of the K-S statistic D 
in constructing a confidence band on any continuous cumulative distribution 
function centered around the EDF. 


Abernethy et al . (ref. 36) extended the procedure by modeling the sample 
and its parent population with the Weibull distribution. Equation (32), with 
tabular data of K(N) from reference 36, is used in calculating the 90-percent 
confidence bands about the Weibull distribution for complete samples only. 

Thus , 
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Confidence bands = [ F ( x ) - K(N), F<x) + K(N)] 


(32) 


where F(x) is the failure probability obtained by substituting the MLE's of 
the Weibull parameters in equation (4). The Kanofsky functions, denoted by 
K(N), depend on the sample size. He have developed two sample-size-dependent 
approximate polynomial functions to interpolate these values, and we have incor- 
porated them into the SCARE program. The two terms in equation (32) give the 
lower and upper 90-percent confidence band values about the Weibull distribu- 
tion, respectively. Figure 7 indicates that for HIPped SiC flexure bars, the 
fracture stresses lie between the confidence bands about the Weibull line. 

The value of the statistic L is 0.0829. It is also clear that the stress 
data fits a single Weibull distribution quite well. 


Batdorf-Wei bul 1 Statistical Parameters for Reliability Analysis 

After evaluating the parameters C and m (either by the least-squares 
analysis or by the maximum likelihood method) using data from four-point bend 
bars or pure tensile specimens and after assessing the intrinsic uncertainties, 
the parameters required in reliability calculations are determined. The 
Weibull modulus m is directly used by all statistical fracture models and, 
consequently, requires no additional calculations. However, for the Weibull 
model the scale parameter a 0 is needed, while for the Batdorf approach the 
flaw-density coefficient kg must be calculated. Both of these material con- 
stants are functions of m, C, and specimen geometry. Details of evaluating 
a 0 and kg are summarized elsewhere (refs. 11 and 12), but the results are 
included here to show how the two sets of constants are related. For volume- 
flaw analysis, using quadrant (Li = 2 L 2 ), four-point MOR bar data with known 
geometry, a 0 \j is calculated from 


oV 


(L-j + m^l^) 


Cy(m v + 1 )' 


/m. 


(33) 


For shear-insensitive fracture criterion, kgv is obtained from 

,2 


kg V = (2niy + 1 ) 


2C y (m v + 1>‘ 


wh(L 1 + m^l^) 


(34) 


where w is the width of the rectangular beam, h is the height, and L] and 
1_2 are, respectively, the lengths between symmetrically placed outer loads and 
Inner loads. By comparing equations (33) and (34), we concluded that when the 
normal stress failure criterion is used, kgv and a 0 \j are related by 



However, for shear-sensitive fracture criteria when crack shape must also be 
specified, the format of equation (35) changes according to the criterion and 
the crack shape selected. In addition, if we note that C\j = V e a 0 v , where 
V e is the effective volume, equation (33) can be solved for V e to obtain 
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(36) 
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(L^ + m^i^) 
(m v + l ) 2 


(m v + 2) 
4(m y + 1> 2 


where V = whL] , the total volume of the specimen. For uniaxial tensile load- 
ing, the effective volume is equal to the specimen gage volume Vg, which is at 
the same time the total specimen volume V of interest in that test. 

In surface-flaw-based reliability analysis, the effective area A e for a 
rectangular beam in four-point bending is (ref. 12) 


A e = 



' rruw 
wTh + 


1 }) 2L, (w + h) 


(37) 


In terms of A e and the extreme fiber fracture stress of, we can rewrite equa- 
tion (4) as 


P 


fS " 



(38) 


Comparing equations (4) and (38) gives 


a oS ' 



(39) 


Finally, the Batdorf surface-crack density coefficient kgs Is calculated from 
(ref. 12) 


3/2-m, 


‘BS 


c s m sV m s 




(1 + n> 


m S 2 


(1 - n ) 1/2 d n 


-1 


(40) 


where n = [(?of/a cr ) - 1]. Equation (40) can be integrated numerically by 
using Gaussian quadratures to evaluate kgs- In pure tensile loading, the 
effective area in equations (39) and (40) is just equal to the specimen gage 
area Ag. 

Recently, it has been noted that equation (40) can be evaluated in closed 
form to obtai n 
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BS " 
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(41) 
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where r is the gamma function, and A e and Cg were eliminated through the 
use of equation (39). The last term in equation (41) is often referred to as 
the uniaxial Weibull surface-crack coefficient, that is k^g = o 0 S^- To 
derive equations (40) and (41), we have assumed that the material is shear- 
insensitive. Under this condition, the Batdorf flaw-density constant kgs Is 
identical with the polyaxial Weibull coefficient k^pg that can be obtained 
from using the Weibull normal-stress-averaging method. In the previous equa- 
tions, the parameters m and C are the biased MLE's of the parameters them- 
selves for both volume- and surface-flaw fracture characterization. In the 
SCARE computer program, equations (35) and (41) are used to calculate the mate- 
rial constants kgv and kgg, respectively, for the Bartdorf shear-insensitive 
model . 

However, if we select shear-sensitive fracture criteria to convert MOR bar 
or tensile fracture data into Batdorf-model crack-function coefficients, equa- 
tions (35) and (41) will no longer apply. As an example, if we choose to model 
both volume and surface imperfections as Griffith cracks (usually conservative 
but not a good representation of flaws) and apply the total strain energy- 
re 1 ease-rate criterion (self-similar crack extension) to predict impending frac- 
ture, the following equations are obtained: 



In deriving equations (42) and (43), collinear crack growth was assumed, and the 
results for typical values of m y and mg (my = mg = 10.0) give considerably 
smaller Batdorf flaw-density coefficients than those obtained from using the 
previous shear-insensitive theory. Consequently, failure probabilities for a 
given structure will decrease when equations (42) and (43) are used to charac- 
terize a material. Presently, a new effort is underway to derive closed-form 
solutions for the Bardorf coefficients by using out-of-plane crack extension 
criterion (ref. 49) with the penny-shaped crack (PSC) for volume flaws and the 
semicircular-edge crack for surface flaws. In the SCARE program, all these 
user options are included, and results are obtained through numerical proce- 
dures as in multidimensional stress states. It is important to note that using 
the Batdorf theory requires no new material test procedures and data. The 
knowledge of Weibull parameters m and o 0 , whether performing volume- or 
surface-flaw analysis, is adequate for all fast-fracture-failure calculations. 


DISCUSSION AND EXAMPLES 

Several benchmark problems were analyzed from the open literature to calcu- 
late the material strength parameters. The MLE's of the Weibull parameters for 
the results of endurance tests of nearly 5 000 deep-groove ball bearings 
(ref. 37) were compared and are included in table I. The test results, in mil- 
lions of revolutions, for 23 ball bearings were 17.88, 28.92, 33.00, 41.52, 
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42.12, 45.60, 48.48, 51.84, 51.96, 54.12, 55.56, 67.80, 68.64, 68.64, 68 . 88 , 

84.12, 93.12, 98.64, 105.12, 105.84, 127.92, 128.04, 173.40. Table I shows 
that the SCARE results are very close to that of Thoman et al . (ref. 20) and 
Lieblein and Zelen (ref. 37). For the preceding data set, the data value 
173.40 is detected as an outlier at 1-percent significance level, the K-S sta- 
tistic D is 0.1512 (fig. 6 ) at 66.7-percent probability level, and the A-D 
statistic A 2 is 0.3288 at 91.5-percent probability level. This indicates 
that the data fits the hypothesized Heibull distribution quite well, especially 
since only 23 specimens are available in this sample. 

The IEA Annex II agreement (ref. 38) focuses on cooperative research and 
development among several countries in the areas of structural ceramics. In 
November 1986, 400 HIPped SiC flexure bars from Germany were distributed by 
ORNL to the five U.S. participating laboratories, including NASA Lewis Research 
Center. The bars were fractured at these laboratories, and the fracture stress 
data sets were returned to ORNL as complete data without censoring for differ- 
ent failure modes. Details of the statistical analyses of these data sets are 
given in reference 38. The results of 80 SiC flexure bars tested at NASA Lewis 
were obtained, (table II) and analyzed by using our present enhancement to 
SCARE to calculate the MLE's of the Heibull parameters. According to the 
results in table III, the Heibull parameter values from SCARE match two other 
predictions very well and are in reasonable agreement with the other three pre- 
dictions, despite the use of different test operators, test equipment, and soft- 
ware to reduce the data. The only difference between the SCARE and NASA Lewis 
data shown is the use of different software to calculate m and ctq. ASEA 
Cerama HIPped S i 3 N 4 bars (ref. 38) from Sweden were also fractured at NASA 
Lewis, and, subsequently, the shape parameter value was independently obtained 
by using the least-squares analysis and the maximum likelihood methods. This 
data set has also been analyzed by using the SCARE estimators, and the compari- 
son of results is shown in table IV. Finally, several proprietary Monte Carlo 
simulated distributions were analyzed, and the results were consistently in 
good agreement with those from other estimators. 


CONCLUSIONS 

The general purpose, statistical, fast-fracture probability code SCARE has 
been enhanced to include the details of calculating the shape and scale parame- 
ters of the two-parameter Heibull distribution by using the least-squares analy- 
sis and the maximum likelihood methods for volume- and surface-flaw-induced 
fracture in ceramics. If any sample, complete or censored, fits the Heibull 
distribution, hhen these two methods are the most powerful and appropriate tech- 
niques for obtaining the Heibull parameter estimates. In addition, for com- 
plete samples only, the maximum likelihood method enables the user to evaluate 
90-percent confidence intervals for the MLE's of the Heibull parameters and the 
unbiased estimates of the shape parameter. Therefore, further research is 
needed in developing confidence intervals and the unbiasing factors for cen- 
sored samples. The current addition to SCARE describes how to detect outliers 
and how to calculate Kolmogorov-Smirnov and Anderson-Darling goodness-of-f i t 
test statistics and the 90-percent confidence bands about the Heibull 
distribution . 
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The fracture data of a material, especially in the limited-strength range 
used in laboratory tests, may fit many distributions other than the Weibull dis- 
tribution. However, prior experience with ceramics suggests that the Weibull 
distribution is the best distribution to fit the failure data. Fractography is 
essential for proper identification of flaw origins. Independent Weibull equa- 
tions are used for volume- and surface-flaw-based reliability analyses. In 
addition, the Weibull parameters are used to calculate the Batdorf flaw-density 
constants for shear-insensitive cracks using closed-form solutions. These con- 
stants are calculated for both volume- and surface-flaw-based reliability analy- 
sis. For shear-sensitive cracks, however, numerical procedures are used in the 
SCARE code, just as in multidimensional stress states. 
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TABLE I . - ENDURANCE TESTS OF DEEP-GROOVE 


BALL BEARINGS 

[All estimates are biased estimates.] 


Parameter 

estimates 

Maximum 1 ikel ihood 
estimates (MLE) of - 


Weibul 1 
modulus , 
m 

Characteri stic 
strength, 

0e ’ K 
Revolution x 10~° 

Lieblein and Zelen 
(ref. 37) 

2.230 

80.00 

Thoman et al . 
(ref. 20) 

2.102 

81 .99 

SCARE 

(refs. 11 and 12) 

2.103 

CO 

CO 

CO 
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TABLE II. - EXTREME FIBER FRACTURE STRESSES OF ELECKTROSCHMELZWERK KEMPTEN (ESK) 

HIPped SILICON CARBIDE (SIC) BARS 


Flexure 

bar 


Strength , 
MPa 


41 

516.2 

42 

519.8 

43 

527.6 

44 

530.7 

45 

530.7 

46 

545.7 

47 

548.8 

48 

552.7 

49 

559.6 

50 

562.4 

51 

563.3 

52 

566.1 

53 

566.5 

54 

570.1 

55 

572.8 

56 

575.0 

57 

576.1 

58 

580.0 

59 

582.6 

60 

588.0 


TABLE III. - WEI BULL PARAMETERS FOR FLEXURE 
STRENGTH DATA SET OF ESK HIPped 
SILICON CARBIDE BARS 
[All estimates are biased estimates.] 


Fracture data 
set 


Maximum like! ihood 
estimates (MLE) of - 


Wei bul 1 
modulus , 
m 


Characteri Stic 
strength, 


541 .17 

517.01 
508.88 
556.72 

554.01 
555.80 






TABLE IV. - WEIBULL PARAMETERS FOR FLEXURE STRENGTH DATA SET 
OF ASEA CERAMA HIPped SILICON NITRIDE BARS 
[All estimates are biased estimates.] 



FIGURE 1. - FLOWCHART FOR MATERIAL STATISTICAL STRENGTH PARAMETERS CALCULATIONS. 
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PROBABILITY DENSITY FUNCTION 


ORIGINAL PAGE IS 
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FIGURE 2. - PROBABILITY DENSITY FUNCTION- STRENGTH DIAGRAM 
FOR VARIOUS VALUES OF THE WEIBULL SHAPE PARAMETER. 



NUMBER OF SAMPLES 

FIGURE 4. - WEIBULL SLOPE ERROR VERSUS SAMPLE SIZE. 



FRACTURE STRESS, MPa 

FIGURE 3. - COMPLETE SAMPLE AND ASSUMED CENSORED 
SAMPLE WEIBULL DISTRIBUTIONS, AND OUTLIER IN 
ASEA CERAMA HIPped SILICON NITRIDE (Si 3 N^). 
(FRACTURE STRESS DATA GENERATED AT NASA-LEWIS; 
NOT ALL DATA POINTS SHOWN.) 



NUMBER OF SAMPLES 

FIGURE 5. - CHARACTERISTIC VALUE ERROR VERSUS SAMPLE SIZE - 
90% CONFIDENCE BANDS. 
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PROBABILITY OF FAILURE 



FIGURE 8. - EMPIRICAL DISTRIBUTION FUNCTION CEDF) OF DEEP-GROOVE BALL-BEARING ENDURANCE 
TESTS. (K-S STATISTIC D = D + = 0.1512.) 



FIGURE 7. - 90-PERCENT CONFIDENCE BANDS ABOUT THE 
WEIBULL LINE FOR ESK HIPped SILICON CARBIDE 
(SiC). (FRACTURE STRESS DATA GENERATED AT NASA- 
LEWIS; NOT ALL DATA POINTS SHOWN; L = 0.0829.) 
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